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ABSTRACT 

The properties of uniformly rotating white dwarfs (RWDs) are analyzed within the framework of general 
relativity. Hartle's formalism is applied to construct the internal and external solutions to the Einstein equations. 
The WD matter is described by the relativistic Feynman-Metropolis-Teller equation of state which generalizes 
the Salpeter's one by taking into account the finite size of the nuclei, the Coulomb interactions as well as 
electroweak equilibrium in a self-consistent relativistic fashion. The mass M, radius R, angular momentum /, 
eccentricity e, and quadrupole moment Q of RWDs are calculated as a function of the central density p c and 
rotation angular velocity Q. We construct the region of stability of RWDs (J-M plane) taking into account 
the mass-shedding limit, inverse /3-decay instability, and the boundary established by the turning-points of 
constant J sequences which separates stable from secularly unstable configurations. We found the minimum 
rotation periods ~ 0.3, 0.5, 0.7 and 2.2 seconds and maximum masses ~ 1.500, 1.474, 1.467, 1.202 M for 
4 He, 12 C, ls O, and 56 Fe WDs respectively. By using the turning-point method we found that RWDs can indeed 
be axisymmetrically unstable and we give the range of WD parameters where it occurs. We also construct 
constant rest-mass evolution tracks of RWDs at fixed chemical composition and show that, by loosing angular 
momentum, sub-Chandrasekhar RWDs (mass smaller than maximum static one) can experience both spin-up 
and spin-down epochs depending on their initial mass and rotation period while, super-ChandrasekharRWDs 
(mass larger than maximum static one), only spin-up. 

Subject headings: Hartle's Formalism — Rotating White Dwarfs — Rotational Instability 



1. INTRODUCTION 

The relevance of rotation in enhancing the maximum sta- 
ble mass of a white dwarf (WD) have been discussed for 
many years both for uniform rotation (se e e.g. Jamesl [l964[ 
Anandl 119651: iRoxburgh & Durnevl 119661: [Monaghan 119661: 
Geroyannis & H adionoulos 1989) and differential rotation 
(see e.g. lOstriker & B odenheimer 1968; Ostriker & Tassoull 



119691: iTassoul & Ostrikerill970t iDurisenl 11975b . Newtonian 
gravity and post-Newtonian approximation have been mainly 
used to compute t he structure of the star, with the excep- 
tion of the work of lArutyunyan et. aD (119711) . where rotating 
white dwarfs (RWDs) were computed in full General Relativ- 
ity (GR). From the microscopical point of view, the equation 
of state (EOS) of cold WD matter has been assumed to be 
either the one of a microscopically unifo rm degenerate elec- 
tron fluid used by Chandrasek han (119311) in his classic work, 
or assumed to have a polytrop ic form. 

However, as shown first bv ISalpeterl (1196 11) in the Newto- 
nian case and then by (Rotondo et al. 201133) in General Rel- 
ativity (GR), a detailed description of the EOS taking into ac- 
count the effects of the Coulomb interaction are essential for 
the determination of the maximum stable mass of non-rotating 
WDs. Specific microphysics of the ion-electron system form- 
ing a Coulomb lattice, together with the detail computation of 
the inverse yS-decays and the pycnonuclear reaction rates, play 
a fundamental role. 

A new EOS taking into account the finite size of the nu- 
cleus, the Coulomb interactions, and the electroweak equi- 
librium in a self-co nsistent relativistic fas hion has been re- 
cently obtained by iRotondo et al.l (1201 lbl) . This relativistic 
Feynman-Metropolis-Teller (RFMT) EOS generalizes both 



the IChandrasekharl (119311) and ISalpeterl (119611) works in that 
a full treatment of the Coulomb interaction is given through 
the solution of a relativistic Thomas-Fermi model. This 
leads to a more accurate calculation of the energy and pres- 
sure of the Wigner-Seitz cells, hence a more accurate EOS. 
It has been shown how the Salpeter EOS overestimates at 
high densities and underestimates at low densities the elec- 
tron pressure. The application of this new EOS to the struc- 
ture of n on-rotating 4 He, 12 C, lfi O and 56 Fe was recently 
done in (IRotondo et alj 1201 lal) . The new mas s -radiu s re- 
lations generalize the wor ks of IChandrase kharl (119311) and 
lHamada & Salpeten (119611) : smaller maximum masses and a 
larger minimum radii are obtained. Both GR and inverse f3- 
decay can be relevant for the instability of non-rotating WDs 
depending on the nuclear composition, as we can see from 
Table \T\ which summarizes some results of IRotondo et al.l 
d2011aS! 



3 / „ 



4 He 


1.56 x 10 10 


GR 


1.40906 


12 C 


2.12 x 10 10 


GR 


1.38603 


16 


1.94 x 10 10 


inverse jS-decay 


1.38024 


56 Fe 


1.18 x 10' 


inverse /3-decay 


1.10618 



TABLE 1 

Critical density and mass for the gravitational collapse of non-rotatin g 
4 He, 12 C, 16 Q and 56 Fe WDs in GR obtained byIRotondo et al] <201 IaI) . 

BASED ON THE RFMT EOS ROTONDO ET AL. I 201 1b). We INDICATE IN THE THIRD 
COLUMN IF THE CRITICAL DENSITY IS DUE EITHER TO INVERSE /3-DECAY OR TO 
GENERAL RELATIVISTIC EFFECTS. 
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We h ere extend the previous results of IRotondo et al.l 
d2011al) for uniformly RWDs at zero temperat ures obeying 
the RFMT EOS. We use the Hartle's approach (iHarflell 19671) 
to solve the Einstein equations accurately up to second order 
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approximation in the angular velocity of the star. We calcu- 
late the mass M, equatorial R eq and polar R p radii, angular 
momentum J, eccentricity e, and quadrupole moment Q, as 
a function of the central density p c and rotation angular ve- 
locity Q of the WD. We construct also RWD models for the 
Chandrasekhar and Salpeter EOS and compare and contrast 
the differences with the RFMT ones. 

We analyze in detail the stability of RWDs both from the 
microscopic and macroscopic point of view in Sec. [3] Be- 
sides the inverse yS-decay instability, we also study the lim- 
its to the matter density imposed by zero-temperature pyc- 
non uclear fusion reactions using up-to-date the oretical mod- 
els dQi^isitd] HH |Yako2iv]iiD HH). The mass- 
shedding limit as well as the secular axisymmetric instability 
boundary are calculated. 

The general structure and stability boundaries of 4 He, 12 C, 
16 and 56 Fe WDs are discussed in in Sec.|4] From the max- 
imally rotating models (mass-shedding sequence), we calcu- 
late in Sec. [5] the maximum mass of uniformly rotating 4 He, 
12 C, le O and 56 Fe WDs for the Chandrasekhar, Salpeter, and 
RFMT EOS, and compare the results with the existing val- 
ues in the literature. We calculate the minimum(maximum) 
rotation period(frequency) of a RWD for the above nuclear 
compositions, taking into account both inverse jS-decay and 
pycnonuclear restrictions to the density; see Sec. [6] 

We discuss in Sec.|7]the axisymmetric instabilities found in 
this work. A comparison of Newtonian and general relativis- 
tic WDs presented in App. [C] show that this is indeed a gen- 
eral relativistic effect. Furthermore, we estimate in App.lDlthe 
accuracy of the "slow" rotation approximation (power-series 
solutions up to order Q 2 ) for the determination of the maxi- 
mally rotating sequence of WDs. In this line, we calculate the 
rotation to gravitational energy ratio and the deviations from 
spherical symmetry. 

In addition, we construct in Sec.[8]constant rest-mass evo- 
lution tracks of RWDs at fixed chemical composition and 
show that RWDs may experience both spin-up and spin-down 
epochs while loosing angular momentum, depending on their 
initial mass and rotation period. 

Finally, in Sec. [9] we outline some astrophysical implica- 
tions of the results presented in this work, which we summa- 
rized in Sec.fTOl 

2. SPACETIME GEOMETRY AND HARTLE'S FORMALISM 



iHartld (11967b described for the first time the structure of 
rotating objects approximately up to second order terms in 
the angular velocity of the star Q, within GR. In this "slow" 
rotation approximation, the solution of the Einstein equa- 
tions in the exterior vacuum can be written in analytic closed 
form in terms of the total mass M, angular momentum / and 
quadrupole moment Q of the star (see App. [A}. The interior 
metric is constructed by solving numerically a system of ordi- 
nary differential equations for the p erturbation functions (see 
lHartldl 1 967t lHartle & Thornd 19681 for details). 

The spacetime geometry up to order Q 2 , with an appropriate 
choice of co ordinates is, in geometrical units c - G - 1, 
described by (lHartlell 19671) 



ds 2 = [e v(r) [l + 2h (r) + 2h 2 (r)P 2 (cos 0)] - co 2 r 2 sin 2 o) dt 2 
+ 2a>r sin 9dtd<p - e y ' 1+2- 



r - M J =°(r) 



- r 2 [1 + 2k 2 (r)P 2 (cos 0)] (d0 2 + sin" 



dr A 
(1) 



where P2(cos 0) is the Legendre polynomial of second order, 
e vir) wAe m = [l-2M J= °(r)/r] A , andM y=0 (r) are the metric 
functions and mass of the corresponding static (non-rotating) 
solution with the same central density as the rotating one. The 
angular velocity of local inertial frames co(r), proportional to 
Q, as well as the functions ho, h%, mo, ma, k 2 , proportional 
to Q 2 , must be calculated from th e Einstein equations (see 
IHartld 19671; lHartle & Thornd 19681 for details); their analytic 
expressions in the vacuum case can be found in App.lAl 

The parameters M, J and Q, are then obtained for a given 
EOS from the matching procedure between the internal and 
external solutions at the surface of the rotating star. The total 
mass is defined by M - M u0 = M J=0 + 5M, where M J=l) is 
the mass of a static (non-rotating) WD with the same central 
density as M J *°, and SM is the contribution to the mass due 
to rotation. 

3. LIMITS ON THE STABILITY OF ROTATING WHITE DWARFS 
3.1. The mass-shedding limit 

The velocity of particles on the equator of the star cannot 
exceed the Keplerian velocity of "free" particles, computed at 
the same location. In this limit, particles on the star's surface 
keep bound to the star only due to a balance between grav- 
ity and centrifugal forces. The evolution of a star rotating at 
this Keplerian rate is acco mpanied by loss o f mass, becom- 
ing thus unstable (see e.g. Stergioulas 2003, for details). A 
procedure to obtain the maximum possible angular velocity 
of the star before rea ching this limit was developed e.g. by 
iFriedman et al.l (119861) . However, in practice, it is less com- 
plicated to compute the mass-shedding (or Keplerian) angular 
velocity of a rotating star, Q J f?°, by calculating the orbital an- 
gular velocity of a test particle in the external field of the star 
and corotating with it at its equatorial radius, r = R eq . 

For the Hartle-Thorne external solutio n, the Keple r ian an - 
gular velocity ca n be writt en as (see e.g. iTorok et al.l (12008): 
iBini et all (120111) and App. |A2l for details) 

= JG-¥- [l - jF,(R eq ) + j 2 F 2 (R eq ) + qF 3 (R eq )] , 

(2) 

where j = cJ/(GM 2 ) and q = c 4 Q/(G 2 M 3 ) are the dimen- 
sionless angular momentum and qu adrup ole moment, and the 
functions F;(r) are defined in App. IA.2I Thus, the numerical 
value of O^? can be computed by gradually increasing the 
value of the angular velocity of the star, Q., until it reaches the 
value Q J f?° expressed by Eq. (fj). 

It is important to analyze the issue of the accuracy of the 
slow rotation approximations, e.g. accurate up to second or- 
der in the rotation expansion parameter, for the description of 
maximally rotating stars as WDs and neutron stars (NSs). We 
have performed in App. [D] a scrutiny of the actual physical 
request made by the slow rotation regime. Based on this anal- 
ysis, we have checked that the accuracy of the slow rotation 
approximation increases with the density of the WD, and that 
the mass-shedding (Keplerian) sequence of RWDs can be ac- 
curately described by the Q. 2 approximation within an error 
smaller than the one found for rapidly rotating NSs, < 6%. 

3.2. The turning-point criterion and secular axisymmetric 
instability 

In a sequence of increasing central density the mass of 
non-rotating star is limited by the first maximum of the M- 
p c curve, i.e. the turning-point given by the maximum mass, 
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dM/dpc = 0, marks the secular instability point and it coin- 
cides also with the dynamical instability point if the perturba- 
tion obe ys the same EOS as of the equilibrium configuration 
(see e.g. IShapiro & Teukolskyll 19831 for details). The situa- 
tions is, however, much more complicated in the case of rotat- 
ing stars; the determination of axisymmetric dynamical insta- 
bility points implies to find the perturbed solutions with zero 
frequency modes, that is, perturbed configurations whose en- 
ergy (mass) is the same as the unperturbed (equil ibrium) one, 
at second order. Howe ver. [Friedman et al.l (Il9 88 ) formulated, 
based on the works of Sorkinj (119811 119821) . a turning-point 
method to locate the points where secular instability sets in 
for uniformly rotating relativistic stars: along a sequence of 
rotating stars with fixed angular momentum and increasing 
central density, the onset of secular axisymmetric instability 
is given by 

Thus, configurations on the right-side of the maximum 
mass of a /-constant sequence are secularly unstable. Af- 
ter the secular instability sets in, the configuration evolves 
quasi-stationarily until it reaches a point of dynamical in- 
stabili ty where g ravitational collapse should take place (see 
Stergioulas 2003, and references therein). The secular insta- 
bility boundary thus separates stable from unstable stars. It 
is worth stressing here that the turning-point of a constant J 
sequence is a sufficient but not a necessary condition for sec- 
ular instability and therefore it establishes an absolute upper 
bound for the mass (at constant J). We construct the boundary 
given by the turning-points of constant angular momentum 
sequences as given by Eq. (0). The question whether dynam- 
ically unstable RWDs can exist or not on the left-side of the 
turning-point boundary remains an interesting problem and 
deserves fu rther attention in vie w of the very recent results 
obtained by iTakami et alJ (1201 lh for some models of rapidly 
rotating NSs. 

3.3. Inverse /3-decay instability 

It is known that a WD might become unstable against the 
inverse /3-decay process (Z, A) — > (Z - 1, A) through the cap- 
ture of energetic electrons. In order to trigger such a process, 
the electron Fermi energy (with the rest-mass subtracted off) 
must be larger than the mass difference between the initial 
(Z,A) and final (Z - 1,A) nucleus. We denote this thresh- 
old energy as 6^. Usually it is satisfied < £^ and there- 
fore the initial nucleus undergoes two succe ssive decays, i.e . 
(Z,A) -> (Z- 1,A) -> (Z- 2, A) (see e.g. lSabeteri (fl%lh : 
Shapi ro & Teukolskvl (fl983l) '). Some of the possible decay 
channels in WDs with the corresponding known experimental 
threshold energies are listed in Table [2] The electrons in 
the WD may eventually reach the threshold energy to trigger 
a given decay at some critical density fP. . Since the electrons 
are responsible for the internal pressure of the WD, configu- 
rations with p > pP dt become unstable due to the softening 
of the EOS as a result of the electro n capture process (see 
lHarrison et all (fl958b : LSalpeterl (fT96l1) for details). In Tabled 
correspondingly to each threshold energy e^, the critical den- 
sity given by the RFMT EOS is shown; see lRotondo et alJ 
d2011al) for details. 

3.4. Pycnonuclear fusion reactions 



4 He — > 3 H + « — > 4n 


7 v ' 

20.596 


1.39 x 10" 


12 C ^12 B ^12 B£ . 


13.370 


3.97 x 10 10 


16 0-» 16 N^ 16 C 


10.419 


1.94 x 10 10 


56p e _,56 Mn ^56 C r 


3.695 


1.18 x 10 9 



TABLE 2 

Onset for the inverse /3-decay of 4 He, 12 C, 1s O and 56 Fe. The experimental 

VALUES OF THE TH RESHOLD E NERGIES S, HAVE BEEN TAKEN FROM TABLE 1 OF 

Audi et al. 12003); see also ( Wapstra & Bos 1977; Shapiro & Teukolsky 
1983). The corresponding critical density pP . are for the RFMT EOS 



(see Rotondo et al. 201 1a) 



In our WD model, we assume a unique nuclear composition 
(Z, A) throughout the star. We have just seen that inverse /?- 
decay imposes a limit to the density of the WD over which the 
current nuclear composition changes from (Z, A) to (Z - 1 , A). 
There is an additional limit to the nuclear composition of a 
WD. Nuclear reactions proceed with the overcoming of the 
Coulomb barrier by the nuclei in the lattice. In the present 
case of zero temperatures T — 0, the Coulomb barrier can be 
overcome because the zero -point energy of the nuclei (see e.g. 
Shapiro & Teukolsl^[1983|) 



na> B 



Ane 2 Z 2 p 



(4) 



-24 



where e is the fundamental charge and M u = 1.6605 X 10 

g is the atomic mass unit. 

Bas ed o n the pycnonucl e ar rates compu ted by IZeldovicE 
(119581) and lCameronl d 1 9591) . ISalpeterl (11961b estimated that in 
a time of 0.1 Myr, 'H is converted into 4 He at p ~ 5 x 10 4 
g cm 3 , 4 He into 12 C at p ~ 8 x 10 8 g cm 4 , and 12 C into 
Mg at p ~ 6 x 10 9 g cm 4 . The threshold density for the 
pycnonuclear fusion of 16 occurs, for the same reaction time 
0.1 Myr, at p ~ 3 x 10 11 g cm 4 , and for 10 Gyr at ~ 10 11 g 
cm 4 . These densities are much higher that the corresponding 
density for inverse /i-decay of 16 0, p ~ 1.9 x 10 10 gcm -3 (see 
Tabled. The same argument applies to heavier compositions 
e.g. 56 Fe; so pycnonuclear reactions are not important for 
heavier than 12 C in WDs. 

It is important to analyze the case of 4 He WDs in detail. At 
densities p pyc ~ 8 x 10 8 g cm" 3 a 4 He WD should have a mass 
M ~ 1.35M (see e.g. Fig. 3 in lRotondo e"t~aT]|201 lab . How- 
ever, the mass of 4 He WDs is constrained to lower values from 
their previous thermonuclear evolution: a cold star with mass 
> 0.5M G have already burned an appreciable part of its He- 
lium content at earlier stages. T hus, WDs of M > 0.5M o w ith 
4 He cores are very unlikely (see Hamada & Salpet erri961l for 
details). It should be stressed that 4 He WDs with M < 0.5M o 
have central densities p ~ 10 6 g cm 4 dRotondo et al.ll201 lab 
and at such densities pycnonuclear reaction times are longer 
than 10 Gyr, hence unimportant. However, we construct in 
this work 4 He RWDs configurations all the way up to their 
inverse /3-decay limiting density for the sake of completeness, 
keeping in mind that the theoretical 4 He WDs configurations 
with M > O.5M could actually not be present in any astro- 
physical system. 

From the above discussion we conclude that pycnonuclear 
reactions can be relevant only for 12 C WDs. It is important 
to stress here that the reason for which the pycnonuclear re- 
action time, rri , determines the lifetime of a 12 C WD is that 



pyc 



reaction times Trjjr < 10 Gyr are achieved at densities ~ 10 
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Fig. 1. — Pycnonuclear reaction times at zero temperature for C+C fusion 
as a function of the density. 

g cirr 3 , lower than the inverse £1 decay threshold density of 
24 Mg, 24 Mg-> 24 Na-> 24 Ne, p ~ 3.2 x 10 9 g cm- 3 (see e.g. 
Salpeten [T96Tl: iShapiro & Teukolsky||1983l) . Thus, the pyc- 
nonuclear I2 C+ 12 C fusion produces unstable 24 Mg that almost 
instantaneously decay owing to electron captures, and so the 
WD becomes unstable as we discussed in Subsec. [3~Tl 

However, the pycnonuclear reaction rates are not known 
with precision due to theoretical and experimental uncertain- 
ties. lHamada & Salpeterl (11961b had already pointed out in 
their work that the above pycnonuclear density thresholds are 
reliable only within a factor 3 or 4. The uncertainties are re- 
lated to the precise knowledge of the Coulomb tunneling in 
the high density low temperature regime relevant to astro- 
physical systems, e.g. WDs and NSs, as well as with the 
precise structure of the lattice; impurities, crystal imperfec- 
tions, as well as the inhomogeneities of the local electron dis- 
tribution and finite temperature effects, also affect the reac- 
tion rates. The energies for which the so-called astrophysical 
S -factors are known from experiments are larger with respect 
to the energies found in WD and NS crusts, and therefore the 
value of the S -factors have to be obtained theoretically from 
the extrapolation of experimental values using appropriate nu- 
clear models, which at the same time are poorly constrained. 
A detailed comparison between the different theoretical meth- 
ods and approximations used for the co mputation of the pyc- 
nonuclear reaction ra tes can be found in (Gasques et al.T l2005b 
lYakovlev et aTll2006l) . 

The S -factors have been computed in (Gasques et al.ll2005t 
lYakovlev et al.l [2006) using up-to-date nuclear models. Fol- 
lowing these works, we have computed the pycnonuclear re- 
action times for C +C fusion as a function of the density as 
given by Eq. (1B31 >. T^jJ. , which we show in Fig.[T] we refer to 
App.lBlfor details. 

We obtain that for t£+ c = 10 Gyr, p pyc ~ 9.26 x 10 y g 

cnr 3 while, for t£+ c = 0.1 Myr, p pyc ~ 1.59 X 10 10 g cm 4 , 
to be compa red with the value p ~ 6 x 10 9 g cirr 3 estimated 
by Salpeterj dl961l) . In order to compare the threshold den- 
sities for inverse y6-decay and pycnonuclear fusion rates, we 
shall indicate in our mass-density and mass-radius relations 
the above two density values corresponding to these two life- 
times. It is important to stress that the computation of the py- 
cnonuclear reactions rates is subjected to theoret ical and ex- 
perimental u ncertainties (see|Gasques et al.l2005l for details). 
For instance, Hamada & Salpeter (1961) stated that these py- 
cnonuclear critical densities are reliable within a factor 3 or 



4. If three times larger, the above value of p pxc for r^ + c =0.1 
Myr becomes p pyc ~ 4.8 x 10 ll) g crrT 3 , larger than the in- 
verse /3-decay threshold density p^ ~ 3.97 x 10 10 g cirr 3 (see 
Table |2). As we will see in Sec. [7] the turning-point con- 
struction leads to an axisymmetric instability boundary in the 
density range p^ =0 = 2.12 x 10 1() < p < p^ g crrT 3 in a 
specific range of angular velocities. This range of densities is 
particularly close to the above values of p pyc which suggests 
a possible competition between different instabilities at high 
densities. 

4. WD STRUCTURE AND STABILITY BOUNDARIES 

The structure of uniformly RW Ds ha v e bee n stud- 
ied by several authors (se e e.g. iJamesI 119641 lAnandl 
1965b Roxburgh & Durneyl 119661: iMonaghanl 119661 
Geroyannis & H adjopoulos 1989). The issue of the sta- 
bility of both uniformly and differentially rotating WDs has 
been studied as well (see e.g.lOstriker & Bodenheimerll 19681 
Ostriker & Tassoull 119691; iTassoul & Ostrikerlll970i iDurisenl 
19751) . All the above computations were carried out within 
Newtonian gravity or at the post-Newtonian approximation. 
The EOS of cold WD matter has been assumed to be either 
the one of a microscopically uniform degenerate electron 
fluid, which we ref er hereafter as Chandrasekhar EOS 
dChmdrasekharTUB, or assuming a poly tropic EOS. How- 
ever, microscopic screening caused by Coulomb interactions 
as well as the process of inverse /?-decay of the composing 
nuclei cannot be pr operly studied within such EOS (see 
lRotondo"eT al. 2011b a, for details). 

The role of gen eral relativistic effects, shown in 
iRotondo et alj d201 lal) . has been neglected in all the above 
precedent literature. The only exception to this r ule is, up to 
our knowledge, the work of Arut yunyan et. al.l (119711) . who 
investigated uniformly RWDs for the Chandrasekhar EOS 
within GR. T hey use an SI 2 approximation foll owing a method 
developed by Sedrakyan & Chubarvan ( 1968), independently 
of the work of Hartle ( 196 7b. A detailed compar ison of our re- 
sults with the ones of Arutyun yan et. al.1 (1197 ll) can be found 
in App.O 

In Figs.|2]-[3]we show the mass-central density relation and 
the mass-radius relation of general relativistic rotating 12 C 
and O WDs. We explicitly show the boundaries of mass- 
shedding, secular axisymmetric instability, inverse /3-decay, 
and pycnonculear reactions. 

Turning now to the rotation properties, in Fig. |4]we show 
the J-M plane especially focusing on RWDs with masses 
larger than the maximum non-rotating mass, hereafter Super - 
Chandrasekhar WDs (SCWDs). It becomes clear from this 
diagram that SCWDs can be stable only by virtue of their non- 
zero angular momentum: the lower-half of the stability line of 
Fig.lH from J = at M/M^° x all the way up to the value of 
/ at M-j*® ~ 1.06M^> determines the critical(minimum) an- 
gular momentum under which a SCWDs becomes unstable. 
The upper half of the stability line determines, instead, the 
maximum angular momentum that SCWDs can have. 

5. THE MAXIMUM MASS 

The maximum masses of rotating WDs belongs to the Kep- 
lerian sequence (see Figs.|2]-|4]i and it can be expressed as 

M J *° = kM J=0 (5) 

where M J „^\ is the maximum stable mass of non-rotating WDs 
and k is a numerical factor that depends on the chemical com- 
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Fig. 2. — Mass in solar masses versus the central density for C (left panel) and for O (right panel) WDs. The solid curve corresponds to the mass of 
non-rotating WDs, the Keplerian sequence is the red thick dashed curve, the blue thick dotted-dashed curve is the inverse fi instability boundary, and the green 
thick solid curve is the axisymmetric instability boundary. The orange and purple dashed boundaries correspond to the pycnonuclear densities for reaction times 
Tpy C = 10 Gyr and 0.1 Myr, respectively. All rotating stable WDs are in the shaded region. 
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4 He 
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TABLE 3 

Properties of uniformly rotating general relativistic 4 He, 12 C, 16 and 56 Fe WDs: P m j*o is the central density in g cm" 3 corresponding to the rotating 

MAXTMUM MASS M^*®; k IS THE DIMENSI ONLESS FACTOR USED TO EX PRESS THE ROTATING MAXIMUM MASS M^*® AS A FUNCTION OF THE NON-ROTATING MAXIMUM MASS M^° x OF 
WD.S, IN SOLAR MASSES, OBTAINED IN RoTONDO ET AL. I 201 1 A), AS DEFINED IN Eq. (5); THE CORRESPONDING MINIMUM RADIUS IS R M J=o, IN KM; P„„„ IS THE MINIMUM 
ROTATION PERIOD IN SECONDS. We RECALL THAT THE CONFIGURATION WITH P„„„ IS OBTAINED FOR A WD ROTATING AT THE MASS-SHEDDING LIMIT AND WITH CENTRAL DENSITY 
EQUAL TO THE CRITICAL DENSITY FOR INVERSE /J-DECAY (SEE TABLELAND THE RIGHT PANEL OF FlG.|6j. THE POLAR Rp""" AND EQUATORIAL K e T" RADH OF THE CONFIGURATION 
WITH P m j„ ARE ALSO GIVEN IN KM. THE QUANTITY (T/\W\) Pl "'" IS THE RATIO BETWEEN THE KINETIC AND BINDING ENERGIES, THE PARAMETER £ P """ IS THE ECCENTRICITY OF THE 
STAR, ROTATFWG AT P„„„. FrNALLY, f mi " AND q P <»>" ARE THE DIMENSIONLESS ANGULAR MOMENTUM AND QUADRUPOLE MOMENT OF WDs, RESPECTIVELY. 



position, see Table [3]for details. For 4 He, 12 C, ie O, and 56 Fe 
RWDs, we found ~ 1.500, 1.474, 1.467, 1.202 M , re- 
spectively. 

In Table [4] we compare the properties of the configura- 
tion with maximu m mass using different EOS, namely Chan- 
drasekhar yu = 2 (Boshka yev et. alJl201 lL see e.g.), Salpeter, 
and RFMT EOS. A comparison with classical results obtained 
with different treatments and EOS can be found in App.lCl 

It is worth mentioning that the maximum mass of RWDs 
is not associated with a critical maximum density for gravita- 
tional collapse. This is in contrast with the non-rotating case 
where the configuration of maximum mass (turning-point) 



corresponds to a critical maximum density over which the 
WD is unstable against gravitational collapse. 

The angular momentum / along the mass-shedding se- 
quence is not constant and thus the turning-point criterion (01 
does not apply to this sequence. Therefore the configuration 
of maximum rotating mass (0 does not separate stable from 
secular axisymmetrically unstable WDs. We have also veri- 
fied that none of the RWDs belonging to the mass-shedding 
sequence is a turning-point of some J =constant sequence, 
and therefore they are indeed secularly stable. We therefore 
extend the Keplerian sequence all the way up to the critical 
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Nuclear Composition 


EOS 


p M m (g/cm 3 ) 


M J *° 

Rp"" x (km) 


R^" ax (km) 


Mi*°/M 


P M?*° X (sec) 


H = 2 


Chandrasekhar 


1.07 x 10 10 


1198.91 


1583.47 


1.5159 


0.884 


4 He 


Salpeter 
RFMT 


1.07 x 10 10 
5.46 x 10 9 


1193.08 
1458.58 


1575.94 
1932.59 


1.4996 
1.5001 


0.883 
1.199 


12 C 


Salpeter 
RFMT 


1.08 x 10 10 
6.95 x 10 9 


1183.99 
1349.15 


1564.16 
1785.98 


1.4833 
1.4736 


0.878 
1.074 


16 Q 


Salpeter 
RFMT 


1.09 x 10 1U 
7.68 x 10 9 


1178.88 
1308.09 


1556.68 
1730.65 


1.4773 
1.4667 


0.875 
1.027 


56 Fe 


Salpeter 
RFMT 


1.14 x 10 9 
1.18 x 10 9 


2002.43 
2000.11 


2693.17 
2686.06 


1.2050 
1.2017 


2.202 
2.195 



TABLE 4 ^ uo 

The maximum rotating mass of general relativistic uniformly rotating 4 He, 12 C, lfi O and 5s Fe WDs for different EoS. P m j*o , R 1 ^""" , R^j"" , and P M nax ARE 

CENTRAL DENSITY, POLAR AND EQUATORIAL RADII, AND ROTATION PERIOD OF THE CONFIGURATION WITH THE MAXIMUM MASS, M?,*®.. 



density for inverse ft decay, (P^ , see Table [2] and Fig. [2] 

6. THE MINIMUM ROTATION PERIOD 

The minimum rotation period P mi „ of WDs is obtained for 
a configuration rotating at Keplerian angular velocity, at the 
critical inverse /J-decay density; i.e. is the configuration ly- 
ing at the crossing point between the mass-shedding and in- 
verse jS-decay boundaries, see Figs. [2] and |4] For 4 He, 12 C, 
16 0, and 56 Fe RWDs we found the minimum rotation periods 
~ 0.28, 0.50, 0.69 and 2.19 seconds, respectively (see Table|3] 
for details). In Table [4] we compare the properties of the con- 
figuration with minimum rotation period using different EOS, 
namely Chandrasekhar p = 2, Salpeter, and RFMT EOS. 

In the case of 12 C WDs, the minimum period 0.50 sec- 
onds have to be compared with the value obtained assum- 
ing as critical density the threshold for pycnonuclear reac- 
tions. Assuming lifetimes r^| f c = 10 Gyr and 0.1 Myr, cor- 
responding to critical densities p pyc ~ 9.26 x 10 9 g cm -3 
and p pyc ~ 1.59 x 10 10 g cm" 3 , we obtain minimum periods 

P'min ~ 0-95 an d 0.75 seconds, respectively. 

It is interesting to compare and contrast some classical 
results with the ones pres ented in this work. Using post- 
Newtonian approximation, Roxburgh & Durnevl (119661) ana- 
lyzed the problem of dynamical stability of maximally rotat- 
ing RWDs, i.e. WDs rotating at the mass-shedding limit. The 
result was a minimum polar radius of 363 km, assuming the 
Chandrasekhar EOS with p — 2. The Roxburgh critical ra- 
dius is rather small with respect to our minimum polar radii, 
see Table [3] It is clear that such a small radius would lead to 
a configuration with the central density over the limit estab- 



lished by inverse /3-decay: the average density obtained for 
the Roxburgh's critical configuration is ~ 1.47 x 10 10 g/cm 3 , 
assuming the maximum mass 1.48M obtained in the same 
work (see TablelClin App.O. A configuration with this mean 
density will certainly have a central density larger than the 
inverse /J-decay density of 12 C and le O, 3.97 x 10 10 g/cm 3 
and 1.94 x 10 10 g/cm 3 , respectively (see Table |2). The rota- 
tion period of the WD at the point of dynamical instability of 
Roxburgh must be certainly shorter than the minimum values 
presented here. 

The above comparison is in line with the fact that we did not 
find any turning-point that cross the mass-shedding sequence 
(see Figs. [2[{3]l. Presumably, ignoring the limits imposed 
by inverse jS-decay and pycnonuclear reactions, the bound- 
ary determined by the turning-points could cross at some 
higher density the Keplerian sequence. Such a configuration 
should have a central dens ity very similar to the one found by 
iRoxburgh & Durnevl (11966b. 

In the work of lArut yunyan et. al.l (119711) the problem of the 
minimum rotation period of a WD was not considered. How- 
ever, they showed their results for a range of central densities 
covering the range of interest of our analysis. Thus, we have 
interpolated their numerical values of the rotation period of 
WDs in the Keplerian sequence and calculated the precise val- 
ues at the inverse yS-decay threshold for 4 He, I2 C, and I6 that 
have p — 2 and therefore in principle comparable to the Chan- 
drasekhar EOS results with the same mean molecular weight. 
We thus obtained minimum periods ~ 0.31, 0.55, 0.77 sec- 
onds, in agreement with our results (see Table |5}. 

It is important to stress that, although it is possible to com- 



On general relativistic uniformly rotating white dwarfs 7 

Nuclear composition EoS pjL (g/cm 3 ) Rp m '" (km) fffj™ (km) M J p *° /M Q P„>,„ (sec) 

/i = 2 Chandra 1.37 X It)' 1 562.79 734.54 1.4963 0.281 

Salpeter 1.37 X 10" 560.41 731.51 1.4803 0.281 
^He RFMT 1.39x10'' 563.71 735.55 1.4623 0.285 

Salpeter 3.88 X 10 lu 815.98 1070.87 1.4775 0.498 
RFMT 3.97x10'° 816.55 1071.10 1.4618 0.501 

Salpeter 1.89 X 10 lu 1005.62 1324.43 1.4761 0.686 
RFMT 1.94x10'° 1005.03 1323.04 1.4630 0.687 

Salpeter 1.14x10' 2002.43 2693.17 1.2050 2.202 

56 Fe RFMT 1.18x10' 2000.11 2686.06 1.2018 2.195 



TABLE 5 

The minimum rotation period of general relativistic rotating 4 He, 12 C, le O and 56 Fe WDs. p ., is the critical density for inverse B decay. M j *° , R n "" n , and 

R e """ ARE THE MASS, POLAR, AND EQUATORIAL RADII CORRESPONDING TO THE CONFIGURATION WITH MINIMUM ROTATION PERIOD, P m i n . 



pare the results using the Chandrasekhar EOS p = 2 with 
the ones obtained for the RFMT EOS, both qualitative and 
quantitative differences exist between the two treatments. In 
the former a universal mass-density and mass-radius relation 
is obtained assuming p = 2 while, in reality, the configura- 
tions of equilibrium depend on the specific values of Z and 
A in non-trivial way. For instance, 4 He, 12 C, and 16 have 
p — 2 but the configurations of eq uilibrium are rather differ- 
ent. This fact was emphasized by Hamada & Salpeterl (11961b 
in the Newtonian case and further in GR by Rotondoetal. 
d201 lal) . for non-rotating configurations. In Fig.|5]we present 
a comparison of the mass-density and mass-radius for the uni- 
versal Chandrasekhar p = 2 and the RFMT EOS for specific 
nuclear compositions. 

7. OCCURRENCE OF SECULAR AXISYMMETRIC INSTABILITY 

Regarding the s t ability of rotating WDs , 
Ostriker & Bod enheimerl (119681) : lOstriker & Tassoull (11969b : 
Durisen (1975) showed that uniformly rotating Newtonian 
polytropes and WDs described by the uniform degenerate 
electron fluid EOS are axisymmetrically stable at any rotation 
rate. In clear contrast with these results, we have shown here 
that uniformly RWDs can be indeed be secularly axisymmet- 
ric unstable as can be seen from Figs. [2]-|4](green boundary). 
We have constructed in App. [C] Newtonian RWDs for the 
Chandrasekhar EOS and compare the differences with the 
general relativistic counterpart. Apart from the quantitative 
differences for the determination of the mass at high densi- 
ties, it can be seen from Fig. [7] (left panel) the absence of 
turning-points in the Newtonian mass-density relation. This 
can be understood from the fact that the maximum stable 
mass of non-rotating WDs is, in the Newtonian case, reached 
formally at infinite central density. We should then expect 
that turning-points will appear only from a post-Newtonian 
approximation , where the critical mass is s hifted to finite den- 
sities (see e.g. Roxburgh & Durnev 1966, for the calculation 
of dynamical instability for post-Newtonian RWDs obeying 
the Chandrasekhar EOS). 

In this respect the Fig. [4] is of particular astrophysical rel- 
evance. Configurations lying in the filled region are sta- 
ble against mass-shedding, inverse /3-decay and secular ax- 
isymmetric instabilities. RWDs with masses smaller than the 
maximum non-rotating mass (Sub-Chandrasekhar WDs), i.e. 
M J *° < M J n ^\, can have angular momenta ranging from a 
maximum at the mass-shedding limit all the way down to the 
non-rotating limit J — 0. SCWDs, however, are stabilized 
due to rotation and therefore there exist a minimum angular 
momentum, /„„■„ > 0, to guarantee their stability. We have 
shown above that secular axisymmetric instability is relevant 



for the determination of this minimum angular momentum 
of SCWDs (see green boundary in Fig. 0). It is interesting 
to note in this respect that from our results it turns out that 
SCWDs with light chemical compositions such as 4 He and 
12 C, are unstable against axisymmetric, inverse yS-decay and 
mass-shedding instabilities. On the opposite, in SCWDs with 
heavier chemical compositions, such as 16 and 56 Fe, the sec- 
ular axisymmetric instability does not take place; see Fig. [4] 
The existence of the new boundary due to secular axisymmet- 
ric instability is a critical issue for the evolution of SCWDs 
since their lifetime might be reduced depending on their ini- 
tial mass and angular momentum. 

From the quantitative point of view, we have found that ax- 
isymmetric instability sets in for 12 C SCWDs in the range of 
masses M^° y < M < 1.397M Q , for some specific range of ro- 
tation periods > 1.24 seconds. We can express the minimum 
rotation period that a SCWD with a mass M within the above 
mass range can have through the fitting formula 

(M - M J=0 \~°' 61 

p axi = 0.062 — seconds , (6) 

\ M B ) 

where M^. is the maximum mass of general relativistic 
non-rotating 12 C WD s, M w y =° * 1.386M (see Tableland 
iRotondo et alj d201 lal) ). Thus, Eq. © describes the rotation 
periods of the configurations along the green-dotted bound- 
ary in Figs. |2] |2 and [4] Correspondingly, the central density 
along this instability boundary varies from the critical density 
of static 12 C WDs, p^ =0 = 2.12 X 10 10 g cirT 3 (see TableQ), 
up to the inverse /J-decay density, — 3.97 x 10 10 g cm -3 
(see Table[2]). 

It is important to note that at the lower edge of the den- 
sity range for axisymmetric instability, p^f =0 , the timescale 
of C+C pycnonuclear reactions are t^, + c c as 339 yr (see Fig.[T). 
It becomes then of interest to compare this timescale with the 
corresponding one of the secular axisymmetric instability that 
sets in at the same density. 

The growing time of the secular instability is given by the 
dissipation time that can be driven either b y gravitational ra- 
diation or viscosity ( Chandrasekhar! [T970I) . However, grav- 
itational radiation reaction is expected to drive secular in- 
stabilities for systems with rotational to gravitational energy 
ratio T/\W\ ~ 0.14, the bifurcation poi nt between McClau- 
rin spheroids and Jacobi ellipsoids (see Chandrasekhar] [T970i 
for details). Therefore, we expect gravitational radiation 
to become important only for differentially rotating WDs, 
which can attain more mass and more angular momentum 
(Ostriker & Bodenheimer l 1 9681) . In the present case of gen- 
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Fig. 5. — Mass versus central density (left panel) and mass versus equatorial radius (right panel) for general relativistic WDs using the Chandrasekhar and the 
RFMTEOS. 



eral relativistic uniformly RWDs, only the viscosity timescale 
t v is relevant. A rotating star that becomes secularly unstable 
first evolve with a characteristic time r v and eventually reach 
a point of dynamical instability, thus collapsing within a time 

Tj y „ ~ Q^ 1 ~ y/R 3 /GM < 1 s, where R is the radius of the 
star (see e.g. IStergioulasl2003h . 

T he viscosity time scale can be estimated as t,, = R 2 p/j] (see 
e.g. lLindbloml ll987). where p and j] are the density and vis- 
cosity of the star. The viscosity o f a WD assum ing degenerate 
relativistic electrons is given by (lDurisenlll973l) 



/7//^=4.74xl0- 2 ^^p 5/3 



,2x 10 6 



\2/3 

) +1 



(7) 



where Hr(Z) is a slowly varying dimensionless contstant that 
depends on the atomic number Z and the Coulomb to thermal 
energy ratio 

k B T \ 3 2ZMu J 

where kg is the Boltzmann constant and A ^ 2Z has been 
used. 

The expression (0 is valid for values of F smaller than 
the critical value for crystallization T cry . The critical r cry 
is not well constrained b ut its value should be of the order 
of T c n , ~ 100 (see e.g. D urisenlll973t IShapiro & Teu kolskv 
1983). The critical value T cry defines a crystallization tem- 
perature T crv under which the system behaves as a solid. For 
r cry ~ 100, we have T cry * 8 x 10 7 [p/(10 10 gcirr 3 )] 1/3 K, 
for Z = 6. When F approach es r crv the viscosity can increas e 
drastically to values close to (Ivan H orn 1969; Durisen 19731) 



/Z\ 2/3 

r, cry = 4.0 x 10- 2 (-) p 5/6 exp[0.1(F - T cry )] 



(9) 



For instance, we find that at densities p C ' J t and assuming 

' ' cnt & 

a central temperature T > 0.5T cry with T cry as 10 8 K, the vis- 
cous timescale is in the range 10 < t v < 1000 Myr, where the 
upper limit is obtained using Eq. (O and the lower limit with 
Eq. (0. These timescales are longer than the pycnonuclear 
reaction timescale 7^5- = 339 yr at the same density. So, if 
the pycnonuclear reaction rates are accurate, it would imply 
that pycnonuclear reactions are more important to restrict the 
stability of RWDs with respect to the secular instability. How- 
ever, we have to keep in mind that, as discussed in Sec. [ 



the pycnonuclear critical densities are subjected to theoretical 
and experimental uncertainties, which could in principle shift 
them to higher values. For instance, a possible shift of the 
density for pycnonuclear instability with timescales ~ 1 

Myr to higher values p^,l t c > p^jt ~ > wou ld suggest an inter- 
esting competition between secular and pycnonuclear insta- 
bility in the density range p^ ~° < p < p^ . 

8. SPIN-UP AND SPIN-DOWN EVOLUTION 

It is known that at constant rest-mass Mo, entropy S and 
chemical comp osition (Z, A), the spin ev olution of a RWD is 
given by (see Shapiro & Teukolskv 1990, for details) 



Q V dJ Jm ,s,za 



(10) 



where Q = d£l/dt and E = dE/dt, with E the energy of the 
star. 

Thus, if a RWD is loosing energy by some mechanism dur- 
ing its evolution, that is E < 0, the change of the angular 
velocity Q in time depends on the sign of d£l/dJ; RWDs that 
evolve along a track with dill d J > 0, will spin-down (Q < 0) 
and the ones following tracks with dQ/dJ < will spin-up 
(Q>0). 

In Fig. [6] we show, in the left panel, the Q. =constant and 
/ =constant sequences in the mass-central density diagram 
and, in the right panel, contours of constant rest-mass in the 
O - / plane. 

The sign of dil/dJ can be analyzed from the left panel 
plot of Fig. [6] by joining two consecutive J = constant se- 
quences with an horizontal line and taking into account that J 
decreases from left to right and from up to down. The angu- 
lar velocity SI, instead, decreases from right to left and from 
up to down for SCWDs and, for sub-Chandrasekhar WDs, 
from left to right and from up to down. We note that, in the 
SCWDs region SI = constant sequences satisfy dSl/dp c < 
while, in the sub-Chandrasekhar region, both dSl/dp c < and 
dSl/dp c > appear (see minima). SCWDs can only either 
spin-up by angular momentum loss or spin-down by gaining 
angular momentum. In the latter case, the RWD becomes de- 
compressed with time increasing its radius and moment of 
inertia, and then SCWDs following this evolution track will 
end at the mass-shedding limit (see Fig. |6). Some evolution- 
ary tracks of sub-Chandrasekhar WDs and SCWDs are shown 
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with dCl/dJ < 0, spin-up. 



in the right panel of F i g. [6] I t is appropriate to recall here that 
Shapiro & Teukolskvl (I 1 9901) showed that spin-up behavior by 
angular momentum loss occurs for rapidly rotating Newto- 
nian polytropes if the polytropic index is very close to n — 3, 
namely for an adiab atic index T » 4/3. It was shown ex- 
plicitly by Gerovannis & Papasotiriou (2000) that these con- 
ditions are achieved only by Super-Chandrasekhar polytropes. 

Besides the confirmation of the above known result for 
SCWDs in the general relativistic case, we report here the 
presence of minima d£l/dp c = for some sub-Chandrasekhar 
masses (see e.g. the evolution track of the RWD with M = 
1.38M in the right panel of Fig.|6]l which raises the possibil- 
ity that sub-Chandrasekhar WDs can experience, by angular 
momentum loss, not only the intuitively spin-down evolution, 
but also spin-up epochs. 

9. ASTROPHYSICAL IMPLICATIONS 

It is appropriate to analyze the astrophysical consequences 
of the general relativistic RWDs presented in this work. 

Most of the observed magnetic WDs are massive; for in- 
stance REJ 0317-853 wi th M ~ 1.35M an d B ~ (1.7- 
6.6) x 10 s G (see e.g. iBarstow et all [l99l IKulebiet alj 
2010b); PG 1658 +441 with M ~ 1 31M and B ~ 2.3 x 
10 6 G (see e.g. iLiebert et all [l983t iSchmidt et all QUI); 
and PG 10 31+234 with the highest magnetic field ~ 10 9 
G (see e.g. fSchmidt et al.1 [l986t IKulebiet all 120091) . How- 
ever, they are generally found to be slow rotators (see e.g. 
IWickra masinghe & Ferrario 2000) . It is worth men t ioning 
that it has been recently shown bv lGarcfa-Berro et all (120 1 2|) 
that such a magnetic WDs can be indeed the result of the 
merger of double degenerate binaries; the misalignment of 
the final magnetic dipole moment of the newly born RWD 
with the rotation axis of the star depends on the difference of 
the masses of the WD components of the binary. 

The precise computation of the evolution of the rotation pe- 
riod have to account for the actual value at each time of the 
moment of inertia and the equatorial and polar radii of the 
WD. Whether magnetic and gravitational radiation braking 
can explain or not the current relatively long rotation peri- 
ods of some observed magnetic WDs is an important issue 
that deserves the appropriate attention and will be addressed 
elsewhere. 

Magnetic braking of SCWDs has been recently invoked as 



a possible mechanism to explain the delayed time dist ribution 
of type la supernovae (SNe) (see lllkov & Sokerll2012l for de- 
tails): a type la SN explosion is delayed for a time typical of 
the spin-down time scale tb due to magnetic braking, provid- 
ing the result of the merging process of a WD binary system is 
a magnetic SCWD rather than a sub-Chandrasekhar one. The 
characteristic timescale tb of SCWD has been estimated to be 
10 7 < Tg < 10 ll) yr for magnetic fields comprised in the range 
10 6 < B < 10 8 G. A constant moment of inertia ~ 10 49 g cm 2 
and a fixed critical(maximum) rotation angular velocity 
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have been adopted (lllkov & Sokeril2012l) . 

It is important to recall here that, as discussed in Sec. [8] 
SCWDs spin-up by angular momentum loss, and therefore 
the reference to a "spin-down" time scale for them is just his- 
torical. SCWDs then evolve toward the mass-shedding limit, 
which determines in this case the critical angular velocity for 
rotational instability. 

If we express in terms of £2^ =0 (see App. |A.2t , taking 
into account the values of j and q from the numerical inte- 
gration, we find for RWDs that the Keplerian angular velocity 
can be written as 

crQ-jf , (12) 



7*0 



where the coefficient <x varies in the interval [0.78,0.75] in the 
range of central densities [10 5 , 10 11 ] g cm 4 . It is important 
to mention that the above range of cr hold approximately the 
same independently on the chemical composition of the WD. 
However, the actual numerical value of the critical angular 
velocity, Off , is different for different compositions owing 
to the dependence on (Z,A) of mass-radius relation of non- 
rotating WDs. 

Furthermore, as we have shown, the evolution track fol- 
lowed by a SCWD depends strongly on the initial conditions 
of mass and angular momentum as well as on chemical com- 
position, and evolution of the moment of inertia (see Fig. [6] 
and Sec.[8]for details). It is clear that the assumption of fixed 
moment of inertia / ~ 10 49 g cm 2 , leads to a spin-down time 
scale depending only on the magnetic field strength. A de- 
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tailed computation will lead to a strong dependence on the 
mass of the SCWD; resulting in a two-parameter family of 
delayed times tb(M, B). Detailed calculations of the lifetime 
of SCWDs braking-down due to magnetic dipole radiation are 
then needed to shed light on this important matter. Theoretical 
work along these lines is currently in progress and the results 
will be presented in a forthcoming publication. 

Massive fast rotating and highly magnetized WDs have 
been proposed as an alternative scenario of Soft Gamma Ray 
Repeaters (SGRs) an d Anomalous X-ray Pulsars (AXPs); see 
Malhe iro et al.1 (l2012h for details. Within such scenario, the 
range of minimum rotation periods of massive WDs found in 
this work, 0.3 < P m i„ < 2.2 seconds, depending on the nu- 
clear composition (see Table |5j, implies the rotational stabil- 
ity of SGRs and AXPs, which possess observed rotation peri- 
ods 2 < P < 12 seconds. The relatively long minimum period 
of 56 Fe RWDs ~ 2.2 seconds, implies that RWDs describing 
SGRs and AXPs have to be composed of nuclear composi- 
tions lighter than 56 Fe, e.g. I2 C or I6 0. 

10. CONCLUDING REMARKS 

We have calculated the properties of uniformly RWDs 
within the framework of GR using the Hartle formalism and 
our new EOS for cold WD matter bas ed on the relativisti c 
Feynman-Metropolis-Teller treatment dRotondo et al.l201 lbl) . 
which generali zes previous approaches including the EOS of 
Salpeter (1961). A detailed comparison with RWDs described 
by the Chandrasekhar and the Salpeter EOS has been per- 
formed. 

We constructed the region of stability of RWDs taking into 
account the mass-shedding limit, secular axisymmetric insta- 
bility, inverse yS-decay, and pycnonuclear reaction lifetimes. 
The latter have been computed u sing the updated theo retical 
models of iGasques eFall 00051) : lYakovlev et all (120061) . We 
found that the minimum rotation periods for 4 He, 12 C, I6 0, 
and 56 Fe RWDs are ~ 0.3, 0.5, 0.7 and 2.2 seconds, respec- 
tively (see Table |5). For 12 C WDs, the minimum period 0.5 
seconds needs to be compared with the values P''\ = 0.75 



and 0.95 seconds, obtained assuming as critical density the 
threshold for pycnonuclear reactions for lifetimes T^ v t c c =0.1 
Myr and 10 Gyr, respectively. For the same chemical compo- 
sitions, the maximum masses are ~ 1 .500, 1 .474, 1 .467, 1 .202 
M Q (see Table @). These results and additional properties of 
RWDs can be found in Table[3] 

We have presented a new instability boundary of general 
relativistic SCWDs, over which they become axisymmetri- 
cally unstable. We have expressed the range of masses and 
rotation periods where this occurs through a fitting formula 
given by Eq. (0. A comparison with Newtonian RWDs in 
App.lClshow to the conclusion that this new boundary of in- 
stability for uniformly rotating WDs is a general relativistic 
effect. 

We showed that, by loosing angular momentum, sub- 
Chandrasekhar RWDs can experience both spin-up and spin- 
down epochs while, SCWDs, can only spin-up. These results 
are particularly important for the evolution of WDs whose 
masses approach, either from above or from below, the maxi- 
mum non-rotating mass. The knowledge of the actual values 
of the mass, radii, and moment of inertia of massive RWDs 
are relevant for the computation of delay collapse times in the 
models of type la SN explosions. A careful analysis of all the 
possible instability boundaries as the one presented here have 
to be taken into account during the evolution of the WD at 
pre-SN stages. 

We have indicated specific astrophysical systems where the 
results of this work are relevant; for instance the long rotation 
periods of observed massive magnetic WDs; the delayed col- 
lapse of SCWDs as progenitors of type la SNe; and the alter- 
native scenario for SGRs and AXPs based on massive RWDs. 

We would like to thank the anonymous referees for the 
many comments and suggestions that improved the presen- 
tation of our results. J.A.R. is grateful to Enrique Garcfa- 
Berro and Noam Soker for helpful discussions and remarks 
on the properties of magnetic WDs resulting from WD merg- 
ers and on the relevance of this work for the delayed collapse 
of Super-Chandrasekhar WDs. 
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APPENDIX 

THE HARTLE-THORNE SOLUTION AND EQUATORIAL CIRCULAR ORBITS 

The Hartle-Thorne vacuum solution 

The HT metric given by Eq. (HJ can be written in an analytic closed-form in the exterior vacuum case in terms of the total 
mass M, angular momentum /, and quadrupole moment Q of the rotating star. The angular velocity of local inertial frame s u>(r), 
proportional to Q, and the fu nctions ho, Ii2, mo, m2, k.2, proportional to Q 2 , are derived from the Einstein equations (see iHartlel 
1967; Hartle & Thome 1968, for details). Thus, the metric can be then written as 



2M 



1 +2/fciP 2 (cos6») + 2|l - — I ^-(2cos 2 6»- 1) 



1 J 2 



dt 2 + — sin 2 Odtdcp 
r 



where 



and 



J A I M\ 5 (J - J^/M , / r \ 
ki = — t 1 + — +-- t—QU 1 , k 

Mr 3 \ r 8 M 3 V2 \M / 



1 -2\k 



r 4 



P 2 (cos60-2 1 



2M 
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(Al) 



5 Q - J 2 /M 



J 2 5Q- J 2 /M 
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l --> 4m - 1 ) 
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1/2 
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— In 

2 x- 1 



3x 2 -2 



1 



(x 2 - 1) 



3, x+l 

- In r- 

2 x - 1 (x 2 



3x J - 5* 



I) 2 



(A2) 



are the associated Legendre functions of the second kind, with x = r/M - 1, and P 2 (cos 8) = (l/2)(3 cos 2 6 — 1) is the Legendre 
polynomial. The constants M, J and Q the total mass, angular momentum and ma ss quadrupole moment o f the rota ting object, 
respec tivel y. This form of the metric corrects some misprints of the original paper by fHartle & Thomel (ll968) (see also Be rti et. al.l 
(2004) and lBiniet. all d2009l)). The precise numerical values of M, J and Q are calcualted from the matching procedure of the 
exterior and interior metrics at the surface of the star. 

The total mass of a rotating configuration is defined as M — M J *° = M J=0 + 6M, where M y=0 is the mass of non-rotating 
configuration and 6M is the change in mass of the rotating from the non-rotating configuration with the same central density. It 
should be stressed that in the terms involving J 2 and Q the total mass M can be substituted by M J=0 since 6M is already a second 
order term in the angular velocity. 



Angular velocity of equatorial circular orbits 

The four-velocity u of a test particle on a circular orbit in equatorial plane of axisymmetric stationary spacetime can be 
parametrized by the constant angular velocity O with respect to an observer at infinity 

u = T[d t + QS <fl ], (A3) 

where F is a normalization factor which assures that u"u a = 1. From normalization and geodesies conditions we obtain the 
following expressions for T and Q = u^/u' 

r = ±( gtt + 2n gt f + Q 2 s # r 1/2 , g„, r + 2n gt ^ r + n 2 g m = o, (A4) 

hence, £2, the solution of ( lA41 b. is given by 

, . ~gt<P,r ± ^{gt<p,r) 2 ~ gtt,rg<M>,r 

QkorbiO = —T = , (A5) 

" gW,r 
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where (+/-) stands for co-rotating/counter-rotating orbits, ifi and u' are the angular and time components of the four-velocity, 
and a colon stands for partial derivative with respect to the corresponding coordinate. In our case one needs to consider only 
co-rotating orbits (omitting the plus sign in Q, +or t,(r) = Q. or b(r)) to deter mine the mass shedding (Keplerian) angular velocity on 
the surface of the WD. For the Hartle-Thorne external solution Eq. (lAlb we have 

^orb(r) = [l - JFi(r) + j 2 F 2 (r) + qF 3 (r)] , (A6) 

where j = J/M 2 and q - Q/M 3 are the dimensionless angular momentum and quadrupole moment, 

<M\ 3/2 48M 7 - 80M 6 r + 4AfV - 18MV + 40M 3 r 4 + 10MV + 15Mr 6 - 15r 7 



* = 7 > ^ 



\6M 2 r 4 {r - 2M) 



+ F, 



6M 4 - 8M 3 r - 2M 2 r 2 - 3Mr 3 + 3r 4 15(r 3 - 2M 3 ) 
F 3 = ; F, F — In ■ 



16M 2 r(r-2M)/5 ' 32M 3 r-2M' 

The mass shedding limiting angular velocity of a rotating star is the Keplerian angular velocity evaluated at the equator (r = 
R eq ), i-e. 

Of = n orb (r = R eq ). (A7) 

In the static case i.e. when j - hence q — and 6M = we have the well-known Schwarzschild solution and the orbital 
angular velocity for a test particle Q^f on the surface (r = R) of the WD is given by 



Q y=0 



K 



(A8) 



3 



Weak field limit 

Let us estimate the values of j and q recovering physical units with c and G. The dimensionless angular momentum is 

cJ c aMR 2 Q. (£IR\(GM\ 



J -GM 2 ~G M 2 = a {-)(^R) ' (A9) 

where we have used the fact that J = IQ., with / = aMR 2 , and a ~ 0.1 from our numerical integrations. For massive and fast 
rotating WDs we have (SM)/c ~ 1(T 2 and (GM)/(c 2 R) ~ 10 3 , so j ~ 1. 
The dimensionless quadrupole moment q is 



c 4 Q c 4 BMR 2 (GM\- 2 



where we have expressed the mass quadrupole moment Q in terms of mass and radius of the WD, Q = /3MR 2 , where /3 ~ 10~ 2 , 
so we have q ~ 10 4 . 

The large values of j and q might arise some suspicion on the products jF\, j 2 F2 and qF^ as real correction factors in Eq. (IA6) . 
It is easy to check this in the weak field limit M/r <K 1, where the functions Fj can be expanded as a power-series 

so evaluating at r — R 

lQ.R\lGM\ 112 j alQR\lGM\ 2 

so we finally have jF\ ~ 10 -9 ^ 2 , j 2 Fi ~ 10 -9 , and qFj, ~ 10 -2 . We can therefore see that the products are indeed corrections 
factors and, in addition, that effect due to the quadrupolar deformation is larger than the frame-dragging effect. 

PYCNONUCLEAR FUSION REACTION RATES 

The theoretical framework for the determination of the pycnonuclear reaction rates was developed by Sal peter & van Hornl 
(119691) . The number of reactions per unit volume per unit time can be written as 

1 / \ I/3 
^=Z*^(£,)3.90xl0* 5 ^ 4 exp(-2.638/V4)«nii- 3 s- 1 , A= -j , (Bl) 

where S are astrophysical factors in units of Mev barns (1 barn=10~ 24 cm 2 ) that have to be evaluated at the energy E p given by 
Eq. ©. 
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Fig. 7. — Left panel: Mass versus central density of Newtonian and general relativistic WDs for the Chandrasekhar EOS with fi = 2. Both the non-rotating case 
and the Keplerian sequence are shown. We have stopped the density, just for sake of comparison, at the critical density for the onset of inverse jS-decay of 4 He 
p = 1.39 X 10" g cirT 3 . Right panel: Mass versus central density relation f or general relativistic WDs for the Chandrasekhar EOS with /i = 2 for the static and 
the Keplerian sequence in this work and the one of Arutvunvan et. al. 1 1971). 



Treatment/EOS M^J /M q References 



Newtonian/Chandrasekhar /i = 2 


1.474 


Anand (1965) 


Newtonian/Polytrope n = 3 


1.487 


Roxburgh (1965) 


Post-Newtonian/Chandrasekhar fi = 2 


1.482 


Roxburgh & Durnev (1966) 


GR/Chandrasekhar fx = 2 


1.478 


Arutyunyan et. al. ( 1 97 1 ) 



TABLE 6 

Maximum rotating mass of WDs tn literature. 



For the S -factors we adopt the results of Gasq ues~et al.l (|2005) calculated with the NL2 nuclear model parameterization. For 
center of mass energies E > 19.8 MeV, the S -factors can be fitted by 



S(£) = 5.15xl0 16 exp 



3E°308 

-0.428£ - 



1 + e 0.613(8-£) 



MeV barn, (B2) 



which is appropriate for the ranges of the zero-point energies at high densities. For instance, 12 C nuclei at p = 10 1() g cirT 3 have 
a zero-point oscillation energy E p ~ 34 keV. 
All the nuclei (Z, A) at a given density p will fuse in a time T pyc given by 

T = — = P (B3) 

P} Rpyc AM U Rpy C 

where - p/(AM u ) is the ion-density. Gasqu es~et al.l ((2005) estimated that the S -factors (lB2t are uncertain within a factor 
~ 3.5; it is clear from the above equation that for a given lifetime r pyc such uncertainties reflect also in the determination of the 
density threshold. 

COMPARISON WITH THE NEWTONIAN TREATMENT AND OTHER WORKS 

We have constructed so lutions of the Newtonian equilibrium equations for RWDs accurate up to order Q. 2 , following the 
procedure of H artlel (119671) . In Fig. [7](left panel) we compare these Newtonian configurations with general relativistic RWDs for 
the Chandrasekhar EOS with /j — 2. We can see clearly the differences between the two mass-density relations toward the high 
density region, as expected. A most remarkable difference is the existence of axisymmetric instability boundary in the general 
relativistic case, absent in its Newtonian counterpart. 

Up to our knowledge, the only previous work on RWDs within GR is the one of lArutyunyan et. ail (119711). A method to com- 



pute R WDs configurations accurate up to second order in Q wa s deve loped by two of the auth ors ("see lSedrakyan & Chubaryan 
1968, for details), independently of the work of lHartlel (119671) . In (lArutyunyan et. all 1 19711) . RWDs were computed for the 
Chandrasekhar EOS with p. = 2. 

In Fig. [7] (right panel) we show the mass-central density relation obtained with their method with the ones constructed in this 
work for the same EOS. We note here that the results are different even at the level of static configurations, and since the methods 
are based on construction of rotating configurations from seed static ones, those differences extrapolate to the corresponding 
rotating objects. This fact is to be added to the possible additional difference arising from the different way of approaching the 
order Q 2 in the approximation scheme. The differences between the two equilibrium configurations are evident. 

Turning now to the problem of the maximum mass of a RWD, in TablelClwe present the previous results obtained in Newtonian, 
Post-Newtonian approach and GR by several authors. Depending on their method, approach, treatment, theory and numerical 
code the authors showed different results. These maximum mass of RWDs are to be compared with the ones found in this work 
and presented in Table|4]for the Chandrasekhar fi = 2, Salpeter, and RFMT EOS. 



14 



Boshkayev et al. 



ACCURACY OF THE HARTLE'S APPROACH 

In his classic work. iHartld (0*967) described the slow rotation regime by requesting that fractional changes in pressure, energy 
density, and gravitational field due to the rotation of the star are all much smaller with respect to a non-rotating star with the same 
central density. From a dimensional analysis, such a condition implies 

, / c \ 2 GM J=0 

where M J=0 is the mass of the unperturbed configuration and R its radius. The expression on the right is the only multiplicative 
combination of M, R, G, and c, and in the Newtonian limit coincides with the critical Keplerian angular velocity £2^ =0 given by 
Eq. dA8l >. For unperturbed configurations with (GM)/(c 2 R) < 1, the condition (1DU implies Q.R/c <sz 1. Namely, every particle 
must mo ve at non-relativistic velocities if the perturbation to the original geometry have to be small in terms of percentage. 
Eq. ( ID1I ) can be also written as 

Q. « Cl J =° , (D2) 

which is the reason why it is often believed that the slow rotation approximation is not suitable for the description of stars rotating 
at their mass-shedding value. 

Let us discuss this point more carefully. It is clear that the request that the contribution of rotation to pressure, energy density, 
and gravitational field to be small can be summarized in a single expression, Eq. (ID lb . since all of them are quantitatively given by 
the ratio between the rotational and the gravitational energy of the star. The rotational ene rgy i s T ~ M R 2 £r and the gravitational 
energy is \W\ ~ GM 2 /R = (GM/c 2 R)Mc 2 , hence the condition T/\W\ 1 leads to Eq. (EE" or dD2]i. Now we will discuss the 
above condition for realistic values of the rotational and gravitational energy of a rotating star, abandoning the assumption of 
either fiducial or orde r of magnitude calculations. We show below that the actual limiting angular velocity on the right-hand-side 
of the condition ( ID2I ) has to be higher than the Keplerian value. 

We can write the gravitational binding energy of the star as \W\ = jGM 2 /R and the rotational kinetic energy as T = (l/2)/Q 2 = 
(l/2)aMR 2 Q. 2 , where the constants y and a are structure constants that depends on the density and pressure distribution inside 
the star. According to the slow rotation approximation, T/\W\ <K 1, namely 



W\~ jGM 2 IR ~\2y)\W) ~\2y 



n y=0 



« 1, (D3) 



which can be rewritten in analogous form to Eq. ID2l as 



n« J— Oi =u . (D4) 
V a 

Now we check that the ratio of the structural constants is larger than unity. Let us first consider the simplest example of a 
constant density sphere. In this case a = 2/5 and y — 3/5, so ^2y/a « 1.73, and the condition (1D4I > is Q <g 1 .73Qjf °. If we 
consi der now a more realistic density profile, for instance, a polytrope of index n = 3, we have (see e.g. Sha piro & Teukolskvl 
[1981 

5 - n R 2 R 2 23 

where (r 2 ) = 0.1 1303/? 2 . Therefore we have in this case y — 3/2 and a - 0.075, and so Eq. ( ID4b becomes Q. <k 6.32£2^ =0 . This 
is not surprising since T/\W\ — > 0.025 when Q — > fi^ =0 - 

The above analysis has been done assuming spherical symmetry. When deviations from the spherical shape are taken into 
account, the ratio T/\W\ turn to be even smaller than the previous estimates based on spherical polytropes. Since the equatorial 
radius satisfies R eq > R, at mass-shedding we will have Q < fi£ =0 - In fact, in the Roche model the mass-shedding angular 

velocity is Q^- #0 = (2/3) 3 ^ 2 O^ =0 ~ 0.544Q^ =0 , corresponding to a rotational to gravitational energy ratio T/\W\ ~ 0.0074 (see e.g. 
IShapiro & Teukolsky|[T983l) . 

In our RWDs we have obtained that the mass-shedding angular velocity satisfies ~ 0.75£2{f° at any density; see Eq. (fl2l) . 
Accordingly to this, we show in the left panel of Fig. [8] the ratio T/\W\ for RWDs as a function of the central density for the 
Keplerian sequence. For an increasing central density T/\W\ decreases. On the right panel we have plotted the eccentricity versus 
the central density. For increasing central density the eccentricity decreases, so RWDs become less oblate at higher densities. 

Now we turn to evaluate more specifically the deviations from the spherical s ymmetry. The expansion of the radial coordinate 
of a rotating configuration r(R, 9) in powers of the angular velocity is written as (lHartlel ll967 ) 

r = R + £(R, 6) + 0(Q 4 ), (D6) 

where £ is the difference in the radial coordinate, r, between a point located at the polar angle 9 on the surface of constant density 
p{R) in the rotating configuration, and the point located at the same polar angle on the same constant density surface in the 
non-rotating configuration. In the slow rotation regime, the fractional displacement of the surfaces of constant density due to 
the rotation have to be small, namely %(R, 6)/R <K 1, where %(R, 6) = £j (R) + ^ 2 (^)^2(cos 6) and £j (R) and % 2 (R) are function 
of R proportional to Q 2 . On the right panel of Fig. [9] the difference in the radial coordinate over static radius versus the central 
density is shown. Here we see the same tendency as in the case of the eccentricity, that these differences are decreasing with 
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Fig. 8. — Left panel: rotational to gravitational energy ratio versus the central density for maximally rotating RWDs, calculated with the Chandrasekhar EOS 
[i = 1. Right panel: the eccentricity versus the central density for the same sequence of RWDs. 




p c [ g/cm 3 ] p c [ g/cm 3 ] 

Fig. 9. — Left panel: the rotation parameter normalized to the speed of light versus the central density. Right panel: the difference in the radial coordinate 
over the static radius versus the central density. The solid curve corresponds to the difference between equatorial (9 = rr/2) and static radii and the dashed curve 
corresponds to the difference between polar (6 = 0) and static radii. 

an increasing central density. On the left panel the rotation parameter QR/c versus the central density is shown. Here, with an 
increasing central density the rotation parameter increases. Thus, for higher densities the system becomes less oblate, smaller in 
size with a larger rotation parameter i.e. higher angular velocity. 

In order to estimate the accuracy of the slow rotation approximation for RWDs, based on the above results, it is useful to 
compare all the above numbers with the known results for NSs. For instance, we notice that in NSs QM/c ~ 10 , %(R,0)/R ~ 
10~ 2 and£(R,7T/2)/R - 1CT 1 (see e.g. lBerti et. alj2 004), to be compared with the correspo nding values of RWDs shown in Fig. |9j 
CIR/c < 1(T 2 , ^(R,0)/R ~ 1(T 2 and £(R,n/2)/R ~ 10 _1 . iWeber & Glendenningl (f!99l calculate the accuracy o f the Hartle's 
second order approximation and found that the mass of maximally rotating NSs is accurate within an error < 4%; Benh areit al.l 
( 120051) found that the inclusion of third order expansion O 3 improved the mass-shedding limit numerical values in less than 1% 
for NSs obeying different EOS. On the othe r-hand, it is known that the ratio T/\W\ in the case of NSs is as large as ~ 0.1 in the 
Keplerian sequence (see e.g. Tables 1-5 of iBerti & Stergioulasl (120041) ). Since RWDs have T/\W\ and DJi/c smaller than NSs, 
and SR/R = %/R at least of the same order (see left panel of Fig. |8), we expect that the description of the strucure of RWDs up to 
the mass-shedding limit within the Hartle's approach to have at least the same accuracy as in the case of NSs. 



